Particles Sliding on a Fluctuating Surface: Phase Separation and Power Laws 
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We study a system of hard-core particles sliding downwards on a fluctuating one-dimensional 
surface which is characterized by a dynamical exponent z. In numerical simulations, an initially 
random particle density is found to coarsen and obey scaling with a growing length scale ~ f^''^. The 
structure factor deviates from the Porod law in some cases. The steady state is unusual in that the 
density-segregation order parameter shows strong fluctuations. The two-point correlation function 
has a scaling form with a cusp at small argument which we relate to a power law distribution of 
particle cluster sizes. Exact results on a related model of surface depths provides insight into the 
origin of this behaviour. 



PACS numbers; 05.70.Ln, 05.40.-a, 64.75.-fg, 05.70.Jk 



How do density fluctuations evolve in a system of par- 
ticles moving on a fluctuating surface? Can the combina- 
tion of random vibrations and an external force such as 
gravity drive the system towards a state with large-scale 
clustering of particles? Such large-scale clustering driven 
by a fluctuating potential represents an especially inter- 
esting possibility for the behaviour of two coupled sys- 
tems, one of which evolves autonomously but influences 
the dynamics of the other. Semi-autonomous systems are 
currently of interest in diverse contexts, for instance, ad- 
vection of a passive scalar by a fluid [Q, phase ordering 
in rough films Je| , the motion of stuck and flowing grains 
in a sandpile [pL and the threshold of an instability in a 
sedimenting colloidal crystal [Q. 

In this paper, we show that there is an unusual sort 
of phase ordering in a simple model of this sort, namely 
a system of particles sliding downwards under a gravi- 
tational field on a fiuctuating one-dimensional surface. 
The surface evolves through its own dynamics, while the 
motion of particles is guided by local downward slopes; 
since random surface vibrations cause slope changes, they 
constitute a source of nonequilibrium noise for the par- 
ticle system. The mechanism which promotes clustering 
is simple: fluctuations lead particles into potential min- 
ima or valleys, and once together the particles tend to 
stay together as illustrated in Fig. 1. The question is 
whether this tendency towards clustering persists up to 
macroscopic scales. We show below that in fact the par- 
ticle density exhibits coarsening towards a phase-ordered 
state. This state has uncommonly large fluctuations 
which affect its properties in a qualitative way, and make 
it quite different from that in other driven, conserved, 
systems which exhibit coarsening jsj. 

It is useful to state our principal results at the out- 
set. (1) In an infinite system, an initially randomly dis- 
tributed particle density exhibits coarsening with a char- 
acteristic growing length scale C{t) ^ t^^^ where z is the 
dynamical exponent governing fluctuations of the sur- 
face. For some of the models we study, the scaled struc- 
ture factor varies as with a < 1, which rep- 
resents a marked deviation from the Porod law (a — 1) 
for coarsening systems M. Further, a finite system of size 



L reaches a steady state with the following characteris- 
tics: (2) The magnitude of the density-segregation order 
parameter has a nonzero time-averaged value, but shows 
strong fiuctuations which do not decrease as L increases. 
(3) The static two-point correlation function C(r) has a 
cusp at small values of the scaled separation \r/L\. (4) 
The sizes I of particle clusters are distributed according 
to a power law for large Z, up to an L-dependent cutoff. 
These results are established by extensive numerical sim- 
ulations. Further, the properties (3) and (4) are shown 
to be related, using the independent interval approxima- 
tion applied to the cluster distribution. Also, we de- 
fine a related coarse-grained depth model of the surface, 
and show analytically that the steady state characteris- 
tics (2)- (4) hold for this model. 




FIG. 1. Depicting clustering of particles (•) in a section of 
the fluctuating surface. A surface fluctuation (a) — » (6) causes 
the particles to roll into a valley. Particles remain clustered 
even after a reverse surface fluctuation (b) — > (c) occurs. 

The Sliding Particle (SP) model is defined as a lattice 
model of particles moving on a fiuctuating surface. Both 
the particles and the surface degrees of freedom are rep- 
resented by ±l-valued Ising variables {ci} and {rj_ i } on 
a one-dimensional lattice with periodic boundary condi- 
tions, where a spins occupy lattice sites, and r spins the 
links between sites. Then = ^(1 -I- a.i) represents the 
particle occupation of site i, whereas Tj_i = +1 or — 1 
represents the local slope of the surface (denoted / or \, 
respectively). The dynamics of the interface is that of the 
single-step model M, with stochastic corner flips involv- 
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ing exchange of adjacent r's; thus, /\ — > \/ with ratepi, 
while \/ — > /\ with rate qi. A particle and a hole on adja- 
cent sites + exchange with rates that depend on the 
intervening local slope t^-i; thus the moves •\o — > o\m 
and o/» — !■ •/o occur at rate p2, while the inverse moves 
occur with rate (72 7^ P2- The asymmetry of the rates re- 
flects the fact that it is easier to move downwards along 
the gravitational field. Note that the dynamics conserves 
^ (T and ^ r; we work in the sector where both vanish. 
In the remainder of the paper, we report results for sym- 
metric surface fluctuations {pi = qi), whose behavior at 
large length and time scale is described by the continuum 
Edwards- Wilkinson model Further we consider the 
strong-fleld (52 — 0) limit for the particle system, and set 
P2 = Pi- We have also investigated the behaviour away 
from these limits, and found that our broad conclusions 
remain unaffected. The SP model is a limiting case of 
the Lahiri-Ramaswamy (LR) model of scdimenting col- 
loidal crystals ; it corresponds to the tilt field evolving 
autonomously. As we shall see below, this change causes 
strong phase separation 0| to disappear, and a new type 
of phase ordering, characterized by strong fluctuations, 
to appear in its stead. 

In the SP model, particles preferentially occupy the 
lower portions or large valleys of the fluctuating sur- 
face. In order to study the dynamics of the hills and 
valleys of the surface, we define a height profile {hi} 
with hi ~ X]i<j<i'^j-i- We then define a Coarse-grained 
Depth (CD) model by considering spins Si — —sgn^hi) 
where Si is —1 or if the surface height hi at site 
i is below, above or at the zero level. A stretch of like 
Si's = represents a valley with respect to the zero 
level. The time evolution of the CD model variables {si} 
is induced by the underlying dynamics of the bond vari- 
ables {Tj_i }. The model is similar to the domain growth 
model of Kim et al . 

We studied the evolution of the density in the SP model 
starting from an initial random placement of particles 
on the fluctuating surface. After an initial quick down- 
ward slide into local valleys, the density distribution is 
guided by the evolution of the surface profile. To quan- 
tify the tendency towards clustering, we monitored the 
equal time correlation function C{r,t) = {ai{t)ai+r{t)) 
by numerical simulation (Fig. 2). If 2; is the dynamical 
exponent characteristic of the surface fluctuations (z = 2 
for the symmetric surface model), we expect the scale 
C{t) for density fluctuations to be set by the base lengths 
of typical coarse-grained hills which have overturned in 
time t, i.e. C{t) ~ t^/^. This is indeed the case in the 
scaling limit (r >> 1, t » 1, r/C{t) fixed) as shown by 
the collapse to a scaling function Cs{y = r/C{t)) in Fig. 
2. We have also checked that similar scaling collapses 
occur when we use other models of surface fiuctuations 
with widely different values of z {z ~ 3/2 for Kardar- 
Parisi-Zhang (KPZ) surfaces, and z ~ 4 for the Das 
Sarma-Tamborenea model p^). 

The existence of a single growing length scale £{t) is 



indicative of coarsening towards a phase ordered state 
In the SP model, coarsening is driven by surface 
fluctuations, rather than more customary temperature 
quenches. This causes an interesting feature to appear 
in the scaling function Cs, namely a distinctive cusp for 
small argument: Cs{y) — Cq — Ciy" for y « 1 (Fig. 2). 
We find the cusp exponent a ~ 0.5. This cusp implies 
that the scaled structure factor S ^ (fc£)~(^+") for large 
kC. This is substantially different from the Porod law 
behavior (fc>C)^^, characteristic of customary coarsening 
systems 0. 
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FIG. 2. The data for C{r,t) at different times t = 400x2" 
(with n — 1,...,6) shown in the inset, is shown to collapse 
when scaled by £,{t) ~ t^^'^. 

In the steady state, in order to characterize phase or- 
dering of the SP model, we monitored the magnitude of 
the Fourier components of the density profile 
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|l3i we used the lowest 



Q{^) as a measure 



where to = 1 , . . . , L — 1 . As in 
nonzero Fourier component Q* _ ^yj^ 
of the phase separation in our system with conserved dy- 
namics. The time average (Q*) ~ in a disordered state, 
and is ~ 0.318 in a fully phase separated state. For the 
SP model. Fig. 3 shows numerical results for {Q{k)) ver- 
sus k for various system sizes. While (Q*) approaches 
a finite limit as L — s- 00, the values of Fourier compo- 
nents at fixed k decrease with increasing L. This pro- 
vides strong evidence for phase separation, corresponding 
to the occurrence of density inhomogeneities of the order 
of the system size. However Q*{t) shows strong fluctua- 
tions as a function of time (Fig. 3, inset). With increas- 
ing L, the separations between fluctuations increase with 
L, but the amplitude of fluctuations does not decrease. 
The best way to characterize these strong fluctuations is 
through the full probability distribution Prob{Q*) which 
we found numerically. The mean value (Q*) ~ 0.18, 
while the RMS fluctuation {{Q*^) - {Q*fy/^ ~ 0.07. 
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Despite these large fluctuations in Q* , the configuration 
of the system does not become randomly disordered even 
when Q*{t) is small. Rather, we found that the next few 
Fourier modes (fc = 47r/i, Gtt/L, ...) are excited at such 
times, indicating a break-up into a few macroscopic re- 
gions of different density. 
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FIG. 4. C(r) in the steady state of the SP and the CD 
models for L — 64, 128, 256, and 512. Notice the cusp at 
small \r\/L. The curves are fits to the form Co — ciy^^^ + C2j/. 
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FIG. 3. {Q{k)) as a function of k, for different sizes 
L = 32,64, 128 and 256 (from above to below). Inset: Time 
variation of Q* in the steady state for L — 128. 



Not surprisingly, these macroscopic fluctuations leave 
a strong imprint on spatial correlation functions and 
cluster distributions. For instance, the two-point cor- 
relation function C(r) = (ajaj-i-r) varies with r on the 
scale of the system size L, as is evident from Fig. 4 
which shows that the data collapse onto a single curve 
when plotted versus r/L. For comparison recall that 
a phase-separated system with sharp interfaces between 
two macroscopic phases would show a linear decrease 
C(r) = (1 — 2\r\/L) on length scales larger than the 
correlation length. By contrast, the curve in Fig. 4 is 
nonlinear in the full range of r/L. Further as with the 
coarsening correlation function, the scaling function Cs 
for steady state correlations has a cusp for small argu- 
ment: 
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FIG. 5. Probability distribution P{1) of particle cluster 
lengths for the SP and CD models for L = 256, 512, 1024 
and 2048. The straight lines have slopes —1.8 and —1.5, re- 
spectively. 

The relationship between C(r) and P{1) may be un- 
derstood within the independent interval approximation 
(IIA) [Q. Within this scheme, the joint probability of 
having n successive intervals is treated as the product of 
the distribution of single intervals. In our case, the inter- 
vals are successive clusters of particles and holes, which 
occur with probability P{1). Defining the Laplace trans- 
form P{s) = J^dle~^T{l), and C{s) analogously, we 
have M 



with a ~ 0.5. We will see below that this is related to 
the form of the size distribution of clusters, defined as 
groups of contiguous lattice sites occupied by particles. 
Figure 5 shows that this distribution follows a power law 
P{1) - l-^ with e ~ 1.8. 



s(l-sC(s)) 



2 l-P(s) 
(0 1 + 



(3) 



where (l) is the mean cluster size. In usual applica- 
tions of the IIA, the interval distribution P{1) has a 
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finite first moment (l) independent of L. But that is 
not the case here, as P{1) decays as a slow power law 
P{1) ~ l^^ for I >> 1. Since the largest possible value 
of I is L, we have (1) « aL^~^ for large enough L. Con- 
sidering s in the range 1/L << s « 1, we may ex- 
pand P{s) w 1 — bs^~^; then to leading order, the right 
hand side of Eq. (||) becomes bs^~^ /aL'^~^ , implying 
C{s) « 1/s - b/{aL^~^s^-^). This leads to 



C{r) 



1 - 



ar(3 



(4) 



This has the same scaling form as Eq. (|^) . Matching the 
cusp singularity in Eqs. (H) and (0), we get 



a = 2 



(IIA). 



(5) 



By comparison, the numerically determined values for 
the SP model yield 6* -I- a ~ 2.3. We conclude that the 
IIA provides a useful insight into the relationship between 
cluster statistics and correlation functions, even though 
it is not exact. 

We can gain considerable understanding into the na- 
ture of phase ordering in the SP model by analyzing the 
closely related CD model of surface fluctuations. The re- 
lationship between these models, which is plausible at a 
qualitative level, can be quantified by checking whether 
the overlap O — {(TiSi) is nonzero; numerically we find 
O ~ 0.26. The similarity in the behaviour of the two 
models is brought out in Figs. 4 and 5; C{r/L) has a 
cusp, and P{1) a power law decay, for both models. 

We now show that the CD model can be mapped onto 
a random walk (RW) problem and that the mapping can 
be exploited to give exact results for several properties of 
the model. We make a correspondence between each sur- 
face configurations and an RW trajectory, by interpreting 
'ri-i/2 = +1 or —1 as the rightward or leftward RW step 
at the i'th time instant. Then = 1, — 1 or depending 
on whether the walker is to the right, to the left or at 
the origin after the i'th step. Evidently, the lengths of 
clusters of s = 1 spins (or s — —1 spins) represent times 
between successive returns to the origin. Thus, P{1) for 
the CD model is just the well-known distribution for RW 
return times to the origin, which behaves as Z^'^/^ for 
large /. Thus 6 = 3/2 in this model. 

Since successive RW returns to the origin are indepen- 
dent events, the IIA is exact for the CD model. Thus Eq. 
(H) holds, and we conclude that the correlation-function 
cusp exponent a = 1/2 for this model. 

In systems with strong macroscopic fiuctuations such 
as the SP and CD models, characterization of the phase 
separated steady state requires the full probability dis- 
tribution of the order parameter. For the CD model, 
an appropriate (nonconserved) order parameter is M = 

^ Si , which for the RW represents the excess time a 
walker spends on one side of the origin over the other 
side. In order to respect periodic boundary conditions, 
we need to restrict the ensemble of RWs to those which 



return to the origin after L steps. The full probability 
distribution of M over this ensemble is known from the 
equidistribution theorem on sojourn times of a RW [ p^ : 
Prob{M) = 1/2 for M G [—1,1], i.e. every allowed value 
of M is equally likely. This implies |(M)| = 1/2 and 
{{KP) - (|M|)^)i/2 = 1/VT2. The strong fluctuations in 
M mirror the large fluctuations of Q* in the SP model. 

To summarize, both the SP and CD models exhibit 
a phase ordered steady state with unusual fluctuation 
characteristics, manifested in several related ways: slow 
power law decays of the cluster size distribution, a cusp 
singularity in the scaled two-point correlation function, 
and a probability distribution for the order parameter 
which remains broad even in the thermodynamic limit. 

We have checked numerically that this sort of state 
survives in the SP model even when we depart from the 
strong field limit by allowing (72 to be nonzero, and when 
we vary the ratio P2/P1 of rates of particle and surface 
motion from high to low values. Finally, we found that 
even when we allow different rates pi and qi for upward 
and downward corner flips which makes surface fluctua- 
tions KPZ-like, this type of phase separation persists and 
the scaled correlation function shows a cusp ||l^]. This 
could have interesting implications for the behaviour of 
domain walls in the problem of phase ordering on a rough 
surface The state is destroyed, however, if we allow 
for surface kinematic waves, which transport transverse 
fluctuations of the interface at a finite speed through the 
system. It would be interesting to investigate the long- 
time dynamics of the large-scale fluctuations in the mod- 
els discussed above, and to characterize the steady state 
in higher dimensions. 

We acknowledge very useful discussions with M.R. 
Evans, D. Dhar, S. Ramaswamy, and C. Dasgupta. We 
are especially indebted to S.N. Majumdar for invaluable 
suggestions and criticisms at all stages of this work. 



Rev. Lett. 79, 1150 
Ramaswamy, Phys. 



[1] R.H. Kraichnan, Phys. Rev . Lett. 72, 1016 (19 94). 

[2] B. Drossel and M. Kardar, |cond-mat/0002032 . 

[3] P. Biswas, A. Majumdar, A. Mehta, J.K. Bhattacharjee 

Phys. Rev. E 58, 1266 (1998). 
[4] R. Lahiri and S. Ramaswamy, Phys 

(1997); R. Lahiri, M. Barma and S 

Rev. E 61, 1648 (2000). 
[5] S.J. CorneU and A.J. Bray, Phys. Rev. E 54, 1153 (1996); 

V. Spirin, P.L. Krapivsky and S. Redner, Phys. Rev. E 

60, 2670 (1999). 
[6] A.J. Bray, Adv. Phys. 43, 357 (1994). 
[7] S.N. Majumdar, C. Sire, A.J. Bray, and S.J. CorneU, 

Phys. Rev. Lett. 77, 2867 (1996). 
[8] M. Plischke, Z. Racz and D. Liu, Phys. Rev. B 35, 3485 

(1987). 

[9] S. F. Edwards and D.R. Wilkinson, Proc. R. Soc. Lon- 



4 



don, Scr. A 381, 17 (1982). 
[10] J.M. Kim, A.J. Bray and M.A. Moore, Phys. Rev. A 45, 

8546 (1996). These authors use the instantaneous mean 

height as a reference, while defining s,. 
[11] M. Kardar, G. Parisi and Y.C. Zhang, Phys. Rev. Lett. 

62, 89 (1986). 

[12] S. Das Sarma and P. Tamborenea, Phys. Rev. Lett. 66, 



325 (1991). 

[13] G. Korniss, B. Schrnittniann, and R.K.P. Zia, Europhys. 

Lett. 45, 431 (1999). 
[14] W. Feller, An Introduction to Probability Theory and its 

Applications, (Wiley, New York, 1968), Ghap. Ill, Sec. 9. 
[15] D. Das and M. Barma, unpubUshed. 



5 



